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Abstract 

We report a study of spin dynamics based on simulations with the Polymorphic Tracking 
Code (PTC), exploring the dependence of the static polarization limit on various beam parameters 
and lattice settings for a practical RHIC lattice. 
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INTRODUCTION 


The motion of the spin expectation value (the “spin”) S' of a charged particle traveling in 
the electric and magnetic helds in a circular accelerator is described by the Thomas-BMT 
equation P, dS/dO = VL{z,6) x S, where 6 is the azimuthal angle, z is the location in 
the 6 -dimensional ( 6 D) phase space and O contains the electric and magnetic fields in the 
laboratory frame. For calculations it is convenient to write VL{z,d) = flo{0) +u}{z,6) where 
r 2 o(^) is the contribution from motion on the closed orbit and u}{z,6) is the contribution 
from the synchro-betatron motion. 

It is often necessary to describe spin motion with the help of a unit vector field n{z, 9) 
(the ’’invariant spin held”, or ISF for short) [2] and this object will play a central role in this 
paper. This satishes the Thomas-BMT equation along particle trajectories and is periodic 
with respect to 9\ n{z, 6 + 2n) = n{z, 6). The product Jg = S ■ n is an invariant of motion, 
and the motion of S is simply a precession around the local fi-axis. The spin precession 
frequency around n is characterized by the amplitude dependent spin tune Ug [ 2 ]. Let us 
assume that the 6 D orbital motion is integrable and away from orbital resonances and spin- 
orbit resonances (dehned below). Then the static polarization limit [3] Pum = \{n{z,9))\, 
with the inner average taken over orbital phases, is the maximum achievable equilibrium 
beam polarization on a phase-space torus. On the closed orbit, n is denoted by no and 
it is normally vertical in the arcs and Ug reduces to the closed orbit spin tune z / q . In a 
perfectly aligned planar ring, uq = G 70 , where G = 1.79284739 for protons, and 70 is the 
relativistic factor for the design energy. The ?7-axis diverges from no near the following 
spin-orbit resonances (or “spin resonances” in short), and Pum can be small, 

Ug = k + -t- kyfy + k^n^, k, k^, ky, k^ e Z. (1) 

Note that Ug is usually undehned on orbital resonances, otherwise Ug is only a function of 
the orbital actions J and the optical state of the ring. In particular, if //^(J) is an amplitude 
dependent spin tune, then the fractional part of izUglJ) + lo + I ■ n, with Iq E Tj and I G 
is also an amplitude dependent spin tune. In other words, there is an equivalence class of 
amplitude dependent spin tunes [2]. In addition, the dip in Pum across a spin resonance is 
accompanied by a jump in the amplitude dependent spin tune so that a system can never 
actually sit at a resonance as defined in Eq. (1) |3l Hj. The order of a resonance is defined 
as \kx \ -|- \ky \ + \kz\. Normally Ug stays close to no- 
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In a planar ring, the most important spin resonances are those due to vertical closed 
orbit distortions, which occur near vq = k, /c G Z, namely the imperfection spin resonances; 
and those driven by the vertical betatron oscillations, which occur near z/q = A; ± z/^,, /c G Z, 
namely the hrst order intrinsic spin resonances with \ky \ = 1. The major challenge in a high 
energy polarized proton synchrotron like RHIC is to preserve the beam polarization during 
acceleration mu- The well-known Frossart-Stora formula [7] describes the polarization 
loss after crossing a single isolated spin resonance. Introduction of a pair of diametrically 
opposed orthogonal Siberian snakes [H] renders the closed orbit spin tune to be vq = 0.5 
and independent of the beam energy. Therefore intrinsic resonances are avoided for normal 
Vy during acceleration and even with misalignments, z/q remains close to 0.5 independently 
of the energy so that imperfection resonances are avoided too. However, at rational ver¬ 
tical tunes satisfying the condition 1/2 + k = muy, m,k G Z, there still can be strong 
loss of polarization during acceleration. This phenomeon is traditionally called “snake res¬ 
onance” 010] although for rational Uy the amplitude dependent spin tune does not exist 
so that this condition does not correspond to higher-order resonance as dehned in Eq. . 
Of course in a real ring with misalignments, uq need not be exactly 0.5. Then pairs of real 
resonances in the sense of Eq. [^with irrational Uy can appear, sitting symmetrically on each 
side of the Uy for the “snake resonance”. These doublets can also cause loss of polarization 
during acceleration. 

Inspired by the Froissart-Stora formula, in rings without snakes, it is common practice 
to compute the strengths of the imperfection and intrinsic resonances for a practical lattice, 
and identify the most dangerous ones. Near one of these dangerous spin resonances, one 
can invoke the so-called “single resonance model” mm- This model is based on the 
rotating wave approximation whereby the effect of u is dominated by one particular Fourier 
harmonic ai k = k + kx^x + kyVy -f- kz^zi A:, kx, ky, kz G Z, with \kx\ -l- \ky \ -1- \kz\ < 1, and 
the corresponding resonance strength is e^. This lattice-independent model is analytically 
solvable, and can be extended to include Siberian snakes, often modeled by “point-like” 
spin rotations. These lattice-independent models have been extensively studied for proton 
storage rings. Analytical solutions for z/^ ra and the n-axis [T311T5] have been obtained, as 
well as Piim [16]. Thus as in Ref. na, I', is 0.5 independently of the betatron amplitude, if 
Vy is irrational. Moreover, as shown in naiiH] simulations with these lattice-independent 
models are a great help for understanding the peculiar features of spin dynamics for rational 
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Uy and, in particular, on and near “snake resonances”. 

However, in high-energy proton rings, the basic spin resonances (i.e., without snakes) 
of interest might not be well isolated, then a lattice-dependent study is necessary. For 


example, as we explain in Sections 11A and 11B below, doublets of higher-order resonances 
at irrational Uy can occur near to the Uy of a “snake resonance” since Ug need not be 0.5. 

An extensive lattice-dependent study of the behavior of Pum throughout the whole energy 
range was made in the study of polarized proton beams up to 920 GeV in HERA [3 HH] , 
as an approach complementary to direct spin tracking for acceleration. It is also interesting 
to study Piim for the store conditions of RHIC, with a beam energy of 255GeV, since this 
is relevant for the study of polarization variation during physics stores with constant beam 
energy [2011^. 

In section methods for modeling and simulation based on the Polymorphic Tracking 
Gode (PTG) are described, and in section]^ the results of the simulation are presented. 


I. THE MODEL FOR THE SIMULATIONS AND THE METHODS 

The simulations in this paper utilize the Polymorphic Tracking Gode (PTG) developed 
by E. Forest [22]. Designed to model various geometries of particle accelerators, PTG is 
capable of symplectic tracking of the orbital motion and length-preserving transport of spin 
|23j . where vectors of particle coordinates and Taylor maps can be tracked in a polymorphic 
manner, and where the latter enables the normal form analysis of the one-turn map using 
FPP [211 . Fortran programs have been developed to do the spin tracking, which call PTG 
as a library. The modeling of the RHIG lattice is presented hrst, followed by an explanation 
of the methods of simulation. 

A. Modeling of the RHIC lattice 

The MADX model of the RHIG lattice is exported into an input hie for PTG, which is 
read by the Fortran program. When a particle is tracked through an integration step of 
a magnet body, the orbital transfer map is sandwiched in between two spin kicks in equal 
amounts and each orbital transfer map is a second-order symplectic integrator while each 
spin kick is represented by a 3 x 3 orthogonal rotation matrix. The quadrupoles must be 
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split into many integration steps to ensure the accuracy of orbital and spin tracking. An 
upper limit is set for the spin rotation angle of each integration step, calculated for the 
betatron amplitude of the tracked particle. In this study, we use about 7 integration steps 
for each arc quadrupole, and up to 81 integration steps for the quadrupoles in the hnal focus 
triplets. 

There are two different implementations of a Siberian snake in this study. The hrst 
method implements a zero length spin kick that rotates a spin by 180 degrees around an 
axis in the lattice, namely a “point-like snake”, while the second method implements helical 
dipoles [25] into the lattice, namely a “helical dipole snake”. In PTC, a helical dipole is 
modeled with a symplectic transfer map accurate to the 4-th order. Note that the longitu¬ 
dinal magnetic held inside helical dipoles will introduce a small transverse coupling at large 
orbital excursions. 

This study implements a pair of diametrically opposed snakes and therefore the closed- 
orbit spin tune is 0.5. Note that the spin rotators around IPO (interaction point at 6 o’clock) 
and IP8 (interaction point at 8 o’clock) are not included in this study, but they can also be 
modeled in a similar way. 

The beam-beam interaction is the major beam-current-dependent effect that might affect 
the beam polarization during physics stores. The effect of the beam-beam interaction on 
the beam polarization was studied through long-term tracking, with a lattice-independent 
model |26| and element-by-element tracking [27| , where the beam-beam kick on spin motion 
was taken into account. However, the total beam-beam parameter [2H| from the two IPs 
in RHIC is within -0.015, and the contribution from the linearized beam-beam kick to the 
intrinsic resonance strength, is less than that of an arc quadrupole izg. Therefore, the 
beam-beam spin kick is not too important. Nevertheless, the beam-beam interaction also 
introduces an incoherent tune shift for beam particles. This effect is studied in this work, 
with a thin-lens weak-strong beam-beam kick implemented for the orbital motion. 


B. Simulation method 

The n-axis is calculated in PTC using stroboscopic averaging m- Once n is computed 
for a phase space point z at an azimuth 6, a particle is launched at the same location with 
spin parallel to n and tracked for 5000 turns. If none of the three orbital tunes is rational. 
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then the turn-by-turn orbital coordinates trace phase space points on the same torus, and 
the turn-by-turn spins }, j = 1, 2- ■ -5000 are the local n at the corresponding phase space 
points. Therefore, Pum can be calculated as an average over for such a phase space 
torus. In addition, the amplitude dependent spin tune can be obtained by a Fourier analysis 
of the spin motion according to Ref. [2j. So a particle is launched with its spin perpendicular 
to n, and tracked for 5000 turns, and a turn-by-turn spin series {Sj-} is obtained. The NAFF 
algorithm [30] is then applied to calculate the fundamental frequency of the complex series 
i.e., the amplitude dependent spin tune, in the range [0.4, 0.5]. 


II. SIMULATION OF RHIC POLARIZATION AT STORE 

In this section, the dependence of Pum on various beam parameters is studied for the Run 
12 baseline lattice of the RHIC Blue ring, during physics stores. The default working point 
is = 28.695/29.685, and [3x/[3y = 0.62m/0.64m at the two interaction points IP6 and 
IP8 where the detectors STAR and PHENIX are located, respectively, the chromaticities 
are corrected to ^xHy = 1.85/2.27 with two families of arc sextupoles. n and Pum are 
calculated at IP6, while Vg is independent of the azimuthal angle. The betatron amplitudes 
are normalized with the amplitude corresponding to a normalized 95% emittance of lOvrmm- 
mrad. Two cases with two different beam energies are simulated, and the relevant parameters 
are shown in Table. [1} Note that Case 1 is the beam energy of physics data taking, while 
Case 2 corresponds to the beam energy of a very strong intrinsic resonance around [vy] = 0.7 
during the energy ramp. 


Parameter 

Case 1 

Case 2 

Beam energy (GeV) 

254.8675 

200.6022 

C70 

487 

383.31 

Normalized 95% emittance(7r mm • mrad) 

10.0 

10.0 

Intrinsic resonance strength 

0.0026/0.001P 

0.175 


for G70 = 486.69/487.31, respectively. 

TABLE I; Parameters of the RHIC lattice used in scans of Pum 
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A. Amplitude scan 


For a fixed lattice, different betatron amplitudes contribute to different underlying spin 
resonance strengths. Therefore they will lead to different behaviors of the n-axis on the 
tori and different values of Pum- In this study, a scan of Pum and Ug over different vertical 
betatron amplitudes is investigated for Case 1 and Case 2, while the synchrotron amplitude 
is set to zero. Note that for the cases with helical dipole snakes, the snakes introduce a small 
transverse coupling so that there is a nonzero but small horizontal betatron amplitude. 




(a) Case 1 (b) Case 2 

FIG. 1: (color online) Pum and Vg versus vertical amplitudes for Cases 1 and 2. Here, the 
range of vertical betatron amplitudes is ten times larger in Case 1 than that in Case 2. 

The betatron amplitudes are normalized by IOtt mm ■ mrad. The fractional vertical 
betatron tune near the closed orbit is set to 0.699. The helical dipole snakes are 
implemented for both cases, while the point-like snakes are also implemented for Case 1. 

In addition, the locations of visible higher order spin resonances are indicated as the 
intersections between the horizontal dashed lines and the Ug curves in the plots. Note that 
there are two resonances of the form Vg = Ibvy + integer in Case 1, due to the variation of 

Vy with such large vertical betatron amplitudes. 

In Fig. the behavior of Pum and Vg is compared between Case 1 and Case 2, and between 
different snake implementations. Pum in general becomes smaller with increasing vertical 
betatron amplitude, and has dips near spin resonances. Note that the range of vertical 
betatron amplitudes in Case 1 is 10 times larger than that in Case 2, and Pum decreases 
with amplitude much slower in Case 1 than Case 2, with a much smaller underlying intrinsic 
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resonance strength. In the lattice-independent model with a single vertical resonance driving 
term and two diametrically opposed orthogonal Siberian snakes, it was shown psiin] that i>s 
is 0.5 independently of the betatron amplitude, if the fractional betatron tune is irrational. 
However, when the betatron amplitude becomes larger in the real lattice, the nearby spin 
resonances are no longer isolated in both cases, and this analytical model is violated. Then 
we hnd that h>s is shifted away from 0.5 with amplitude, due to interference between nearby 
spin resonances. Note that the shift of Vg from 0.5 is in general much larger in Case 2 than 
Case 1 for the same betatron amplitude. Moreover, the locations of the spin resonances 
bvy — i's = integer are indicated in the plots and match the sudden dips of Pw^ and the jump 
of Ug. Several other higher order spin resonances are also visible and indicated in the plot, 
and big jumps correlate with wide resonances. In Case 1, the betatron amplitude becomes 
so large that the amplitude dependent orbital tune shift is not negligible, and there are two 
locations corresponding to the same kind of spin resonance Vg = IbUy + integer. In addition, 
in Case 1, except for the location of spin resonances, Pum decreases faster with amplitude 
for the point-like snakes, while Ug shifts faster away from 0.5 with amplitude for the helical 
dipole snakes, and the spin resonance Ug = bUy + integer is wider for the point-like snakes. 
Because the helical dipoles also contribute to the driving term of the spin resonances, for 
this case, it appears that the contribution from the helical dipoles cancels part of the total 
resonance strength mainly driven by the quadrupoles. 


B. Betatron tune scan 

For a hxed vertical betatron amplitude, different vertical betatron tunes correspond to 
different distances from major spin resonances. Therefore they will lead to different values 
of Piim at the same vertical betatron amplitude. In this section, the integer part of betatron 
tunes are kept constant so that “vertical betatron tune” refers to the fractional vertical 
betatron tune. A list of vertical betatron tunes in a selected range is generated with a 
hxed step size, and the lattice is then htted accordingly for each case with a hxed fractional 
horizontal tune = 0.69. is computed for these lattice settings with the same vertical 
betatron amplitude lOvr mm ■ mrad. Helical dipole snakes are used in this simulation. 

Fig. 1^ is a scan in the range Vy G [0.501, 0.98] with a step size of 0.0005. It is clearly seen 
that in Case 2, Pum is generally smaller when Vy is closer to 0.5, which indicates that the 




FIG. 2: (color online) Phm versus fractional vertical betatron tune in the range 
[uy] G [0.501,0.980]. For both cases, the step size of tune scan is 0.0005. The vertical 
normalized emittance is IOtt mm ■ mrad, and the horizontal tune is set to 0.69. Helical 
dipole snakes are implemented in both cases. Two locations of “nonresonant beam 

polarization” are indicated in Case 2 . 


spin resonance = z/y + integer is so strong that it affects the behavior of Pum in the whole 
scan range. In Case 1, however, the strength of the resonance Ug = z/yT integer appears to be 
much smaller. Moreover, several other higher order spin resonances are also visible in Case 2, 
indicating that their widths are comparable or larger than the step size 0.0005 in the vertical 
tune dimension. Several locations of “non-resonant beam depolarization” are also observed 
in this plot, where the the dips in Piim do not correspond to a jump in z/^, i.e., the locations 
of spin resonances. As shown in Ref. [TB] , for the lattice-independent model with an isolated 
vertical resonance driving term and two diametrically opposed orthogonal snakes, Pum can 
be analytically expressed via a special function oq of the resonance strength, betatron tune 
and G 70 , which goes to zero at the locations of “non-resonant beam depolarization” of that 
model. This is an example where the study of these lattice-independent models leads to 
physical insights that are nontrivial to obtain otherwise. 

The vertical betatron tune range [0.67, 0.74] is of particular interests because the vertical 
tunes of current RHIC operations are in this range. The scan result with a step size 0.0001 
is shown in Fig. where helical dipole snakes are applied. In Case 2, with this step size, 
it is clear that the 7/10 “snake resonance” is split into a doublet, due to the fact that Ug 
shifts with Uy. So the locations of the double resonances shift with amplitude as well, while 
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(a) Case 1 (b) Case 2 

FIG. 3: (color online) Pum versus fractional vertical betatron tune in the range of 
[uy] G [0.67, 0.74] for both cases, the step size of the tune scan is 0.0001. Three different 
vertical amplitudes are shown for Case 1 and Case 2. The horizontal tune is set to 0.69. 
Helical dipole snakes are implemented in both cases. 


the widths of these spin resonances increase with amplitude. Due to the transverse coupling 
introduced by the helical dipoles, there is a small blip in P^m at Vy = Vx = 0.69, which 
is not visible if point-like snakes are used instead. However, in Case 1, the widths of the 
spin resonances are very small and invisible in the plot. Moreover, except for the resonance 
locations, it is shown that Pum increases with the vertical tune in this tune range. 


C. Effect of the beam-beam interaction 

A thin-lens weak-strong beam-beam kick to the vertical betatron motion, is implemented 
at IP6, with a vertical beam-beam parameter of —0.012. A beam of 1008 particles, with 
a Gaussian distribution for the vertical coordinates is launched and the Piim is computed 
for each particle’s trajectory (torus). As shown in Fig|^ the calculated turn-by-turn n-axis 
of one particle forms a closed curve on the surface of a unit sphere, and this indicates the 
existence of an n-axis on the particle’s torus in the presence of nonlinear betatron motion. 
Moreover, the effect of beam-beam interaction on Pum is insignificant. Note that the helical 
dipole snakes introduce a small transverse coupling in this example. 
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(a) The n-axis on one phase space torus. (b) Pum and Vy versus vertical amplitude 

FIG. 4: (color online) The effect of beam-beam interaction on Pum- The left plot shows the 
turn-by-turn f?-axis of one tracked particle. The right plot shows P^m versus vertical 
amplitude for the cases with or without beam-beam interaction. In addition, the betatron 
tunes of tracked particles (“Simulation”) are compared with the analytical amplitude 
dependence of the tune shifts (“Analytical”) |2H] , in the presence of the beam-beam 
interaction. This is shown in the small hgure inside the right plot. 1008 particles with a 
Gaussian distribution were launched in the tracking. The beam-beam parameter is —0.012, 

and the fractional tunes are (0.690, 0.695). 


D. Effect of horizontal motion 


The simulations for the cases shown above deal primarily with the vertical betatron 
motion, since the contribution of the small horizontal amplitude introduced by the helical 
dipoles is small for these cases. We now include horizontal motion by initializing a beam 
with a 4D Gaussian distribution. The normalized 95% emittances are lOvr mm ■ mrad for 
both horizontal and vertical planes. Two cases are simulated as shown in Fig.|^ one without 
beam-beam interaction, and the other with a beam-beam kick at IP6, whose beam-beam 
parameter is -0.012. Gompared to the case with only vertical betatron motion, when 4D 
motion is included, Pum spreads out with vertical amplitude to some extent, for different 
trajectories. However, the two cases with or without beam-beam effect do not show much 
difference. This is because the variation of Pum with vertical betatron tune is very small for 
Gase 1 in the tune range between 0.67 and 0.74, as shown in Fig. 
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vertical amplitude(t7y) vertical amplitude(£7y) 

(a) 4D motion without beam-beam interaction (b) 4D motion with beam-beam interaction 

FIG. 5: (color online) Comparison of Pwm versus vertical betatron amplitude, as 
normalized by lOvrmm ■ mrad for 1008 trajectories, between the case with 4D motion, and 
the case with only vertical betatron motion. The left plot shows the case without 
beam-beam interaction, and the right plot shows the case with beam-beam interaction. 
The 95% normalized emittances are IOtt mm ■ mrad for both horizontal and vertical planes. 

A beam-beam kick with a beam-beam parameter -0.012 is implemented at IPO of the 
RHIC lattice, for the right plot, and the tune footprint of the beam particles is shown in a 
small hgure inside the right plot as well. The fractional tunes are set to (0.69,0.695). 

III. CONCLUSION 

In this paper we compute the static polarization limit for a practical RHIC lattice with 
the physics store conditions for various beam parameters, as a step towards understanding 
the polarization evolution at store. All calculations are done on the basis of the Polymorphic 
Tracking Code. It is shown that the vertical betatron oscillation has the dominant effect on 
the behavior of the Pum in contrast with the horizontal oscillations. Note that synchrotron 
motion is not included in this simulation because the synchrotron tune of RHIC is very 
small, namely around 3 x 10““^ at store. In this case the use of stroboscopic averaging to 
hnd the n-axis when synchrotron motion is included requires special studies. Moreover, the 
practical modeling of Siberian snakes with helical dipoles leads to different behavior of Pum 
and Us, in contrast to the implementation of the point-like Siberian snakes. So it is advisable 
to model the snakes carefully. The “nonresonant beam polarization” observed and studied 
in the lattice-independent model is also observed in this lattice-dependent model. Moreover, 
the beam-beam interaction doesn’t have much effect on P\[^ for the parameters under study. 
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In addition, machine imperfections can possibly tilt no from the vertical and thereby lead to 
spin resonances due to horizontal motion. Imperfections can also shift uq away from 0.5 pT] . 
A realistic treatment of various sources of machine imperfections requires a very careful 
lattice modeling, and is beyond the scope of this paper. Nevertheless, this work shows how 
a study of Pum and can give insights at hxed energies that are not available by executing 
simple tracking. 
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